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Abstract 

Transverse momentum broadening of fast partons propagating through a large nucleus 
is proportional to the average color field strength in the nucleus. In this work, the cor- 
responding coefficient is determined in three different frameworks, namely in the color 
dipole approach, in the approach of Baier et al. and in the higher twist factorization 
formalism. This result enables one to use a parametrization of the dipole cross section 
to estimate the values of the gluon transport coefficient and of the higher twist matrix 
element, which is relevant for nuclear broadening. A considerable energy dependence of 
these quantities is found. In addition, numerical calculations are compared to data for 
nuclear broadening of Drell-Yan dileptons, J ftp and T mesons. The scale dependence 
of the strong coupling constant leads to measurable differences between the higher 
twist approach and the other two formalisms. 
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1 Introduction 



A fast parton (quark or gluon) propagating through nuclear matter accumulates transverse 
momentum by multiple interactions with the soft color field of the nucleus. At not too high 
energies, this phenomenon is experimentally accessible by measuring nuclear broadening of 
Drell-Yan (DY) dileptons or of J / ip and T mesons produced in proton-nucleus (jpA) collisions. 
Nuclear broadening is defined as the increase of the mean transverse momentum squared of 
the produced particle compared to proton-proton (pp) collisions, i.e. 



S(Pt) = (Pt) p a- (P 2 T ) 



pp- 



During the past decade, at least three different theoretical approaches have been developed 
to describe this effect, namely the color dipole approach jUE], the approach of Baier et al. 
|Hj and the higher twist factorization formalism j3J E] (see also jE] for earlier work). This 
enormous interest in a QCD based description of nuclear effects is mainly motivated by the 
experimental program at the Relativistic Heavy Ion Collider (RHIC). Data from heavy ion 
collisions at RHIC require a profound theoretical understanding of nuclear effects in terms 
of QCD for a reliable interpretation (see e.g. 

In each of the three approaches [U El El IH Ej , broadening is proportional to a nonper- 
turbative parameter, which has to be determined from experimental data. This often limits 
the predictive power of the theory. Moreover, one would like to know, how the different 
theoretical formulations of transverse momentum broadening relate to one another, and to 
what extend they represent the same (or different) physics. The purpose of this paper is 
to present relations between the nonperturbative parameters, thereby illuminating the con- 
nection between these seemingly very different approaches. Since the nonperturbative input 
to the dipole approach ^ E] is known from processes other than nuclear broadening, one 
can then obtain independent estimates for the parameters of the other two approaches and 
calculate 8(pr) in a parameter free way. In addition, we study the energy dependence of the 
nonperturbative parameters, which has to be known if one wants to extrapolate results from 
fixed target energies to RHIC and the Large Hadron Collider (LHC). 

We shall now briefly summarize the basic formulae for nuclear broadening. In the color 
dipole approach 012], transverse momentum broadening of an energetic parton propagating 
through a large nucleus is given by 

<5<PT)S P oie = 2 PA LC n (0), (2) 

where pA = 0.16 fm -3 is the nuclear density, and L = 3/?a/4 is the average length of the 
nuclear medium traversed by the projectile parton before the hard reaction occurs (Ra is 
the nuclear radius). The index 1Z refers to the color representation of the projectile parton, 
1Z = T for a quark and 1Z = A for a gluon. The nonperturbative physics is parametrized in 
the quantities 

, CU(0) = 9C>(0)/4. (3) 

r T ->-0 



CA0) = ^< (rT) 



T 
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Here, cr^(r T ) is the cross section for scattering a color singlet quark-antiquark (qq) pair with 
transverse separation tt off a nucleon N. This dipole cross section arises from a complicated 
interplay between attenuation and multiple rescattering of the incident parton [2] . 

In the approach of Baier et al. [2j (BDMPS approach hereafter), broadening is related to 
the transport coefficient q-R, 

<*<Pt)bdmps = QnL. (4) 

In this approach, all nonperturbative physics is contained in q-ji, which is a measure for the 
strength of the interaction between the projectile quark and the target. 

The dipole and the BDMPS approach quite obviously describe the same physics, see 
Ref. |H]. Both, Cfc(0) and q-ji can be related to the gluon density of a nucleon. By comparing 
the corresponding expressions in Ref. jH] (for C^(0)) and Ref. [B] (for q-ji), one obtains [TU] . 

q n = 2 PA C 1Z (0)- (5) 

Thus, 5(pt)S P o1c = <*(Pt)bdmps- 

The relation to the higher twist factorization formalism |H |3] is less clear. In the dipole 
and in the BDMPS approach, the projectile parton acquires transverse momentum in a 
random walk through the nuclear medium, thereby undergoing multiple soft rescatterings. 
In the higher twist approach, the (anti-)quark from the projectile proton exchanges only 
one additional soft gluon with the nucleus before the DY dilepton is produced, see Fig. [T] 
Broadening then depends on a particular twist-4 matrix element, which is enhanced by a 
power of A 1 / 3 due to the size of the nucleus (A is the atomic mass of the nucleus) [3], 

= ^s(M*)\l QS A^. (6) 

The quantity Aj^Qg originates from a model of the soft-hard twist-4 matrix element [5], 

TqG(x2) ~ Alqs^ 1 ^ 3 fq/ A ( X2 )' wnere fq/A(x2) is the density of quarks with momentum frac- 
tion X2 in the nucleus. The strong coupling constant a s enters at the characteristic hard scale 
of the process that probes the transverse momentum of the incident parton, i.e. the dilepton 
mass M 2 . Thus, in the higher twist approach, a s is small, even though the exchanged gluon 
is soft. The smallness of a s is crucial for the applicability of the QCD factorization theorem. 




Figure 1: Twist-^ contribution to nuclear broad- 
ening. The projectile antiquark carries momen- 
tum fraction X\ of its parent hadron and under- 
goes one soft rescattering before it annihilates 
with a quark from the nucleus. The Drell-Yan 
process is used to probe the transverse momen- 
tum of the antiquark. 
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2 Relating the dipole approach to the higher twist for- 
malism 



In order to relate all three approaches, one clearly cannot simply set equal the broadening 
given in Eqs. (J2J), (JU) and ©, and then read off a relation between (7^(0), q-ji and Alq S . 
Instead, one has to find a relation between these three quantities in an independent way, and 
only after that, one can check whether all three approaches predict the same (or different) 
5{p\) n . It is then possible to use a model for the dipole cross section to estimate qn and 
Alqs, since is known much better than these two quantities. 
The plan is to relate C^(0) and Al QS to the quantity 

(F 2 ) = ^J dy-(N\F^(y-)F^(0)\N), (7) 

which measures the average color field strength experienced by the projectile parton. In 
Eq. ((7j), P + is the light-cone momentum of the nucleus \A) per nucleon. The index u runs 
over the two transverse directions, and F+ u is the gluon field strength operator. Since we 
are dealing with non-perturbative quantities, the result will of course be model dependent. 

The relation between C^(0) and (F 2 ) can be obtained quite straightforwardly. Note that 
the dipole cross section is related to the gluon density xGn(x) of a nucleon by 

71-2 

(Jgg(x,r T ) = —a s r^xG N (x), (8) 
and in light-cone gauge, the gluon density is given by [TT] . 

xGn{x) = J ^ Q ~' 1XP+V ^ N \ F ^^ F U^)\N). (9) 

What are the relevant scales for xGn{x) and a s l Obviously, the Fourier modes of the 
nuclear color field that give the dominant contribution to broadening are of order 5{p^) n . 
Since this is not much larger than A.q CD in present experiments, the gluon density in this 
case is not a parton distribution like in deeply inelastic scattering (DIS) and should not 
be evolved with the QCD evolution equations. Moreover, the soft gluon carries momentum 
fraction x ~ 5{p\) n / {x\S) of its parent nucleon, which is essentially zero. Here, x\ is the 
momentum fraction of the projectile parton, and S is the hadronic center of mass (cm.) 
energy. We therefore set x — ► in Eq. © and write Cjr(0) as 

CA0) = ^a s (5(p 2 T f)(F 2 ). (10) 

It is important to note that in the dipole approach and in the BDMPS approach, the 
scale of a s is the same as in the gluon density. This is the main difference between these two 
approaches and the higher twist formalism. 
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The next step is to find an expression for A^qg- We shall follow the model assumptions 
of Ref. 0, 



! ) 

— 7 / — 



^.^^ e ix 8 p+v - e(l/ - _ j /l -)e(-yr)^(^l9(o)7 + ?(yo-)^(2/r)^(^)l^) (ii) 

^e^ p+ ^i(A| g -(0) 7 + g ( % -)|^) |^^ e(y --yr)9(-^-)(A|i^(yr)i^ a; (y / r)|A)(12) 
« Xl QS A^f q/A (x 2 ). (13) 

Here, £2 is the momentum fraction of the quark from the nucleus in Fig. ^ The meaning 
of the positions y~ on the light-cone are illustrated in Fig. as well. The step functions 
G ensure that the soft gluon is exchanged before the annihilation. In Eq. (|12jl. the matrix 
element is factorized by introducing an approximate unit operator, 1 « \A)(A\/(2P + V), 
where V is the volume of the nucleus [F2\. I n this step, all correlations between the quark 
and the gluon in Fig. [T]are neglected. As pointed out in Ref. [S], one has \y^\ <C Ra because 
of the rapidly oscillating phase factor in Eq. (fT2|). In addition, \y^ — y'{~\ <C Ra because 
of confinement j^j. This allows one to approximate Q(vq — yf) ~ ©(— yf) ~ 0(— y'^) in 
Eq. ()12jl . With these approximations, the ^-integral factorizes to give the nuclear quark 
density f q /A{x2) ~ Af q /^{x2). The integral over the remaining step function yields a factor 
L, and in the last integration one recovers the right-hand side of Eq. (J7J), though with \N) 
replaced by \A). Assuming that there are no non-trivial nuclear effects on the gluon field, 
the result reads, 

A^ 1/3 = \paL(F>). (14) 

Note that we do not introduce a new model for T g (^2)- Eq- (|14|) follows from the model 
assumptions of Ref. 

Thus, in all three approaches, broadening is related to the quantity (F 2 ), and one finds 
from Eqs. ©, © and ®, 

2vr 2 2tt 2 
<^)£poie = *<p£>£dmps = — a s (5{p 2 T f)p A L(F 2 ) , 5(p 2 T )£ T = —a s (M 2 )p A L(F 2 }. 

(15) 

The new result here is the coefficient 2ir 2 a s pAL/3, the proportionality between broadening 
and the average color field strength in the target was already known before EH E)- It 
is remarkable, that the only difference between S(p 2 -<)^ ipolc and <5(pt)ht * s ^ e sca ^ e °f the 
strong coupling constant. We stress that this difference cannot be dismissed as a higher 
order correction. Instead, it is the result of different physical pictures of nuclear broadening. 

At first sight, the result Eq. (JT5|l may seem puzzling. How can the double scattering 
approximation yield essentially the same expression for broadening as a resummation of all 
rescatterings? In fact, it was demonstrated in Refs. [21 E]' that double scattering does not 
lead to an A -"-^-dependence of broadening. 

This contradiction can be resolved in the following way: The probability to have n 
interactions of the projectile parton with the medium before the Drell-Yan process takes 
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place is (neglecting correlations) Poisson distributed, P n = (aT A ) n e~ aTA /n\, where a is the 
cross section for a single soft scattering, and is the nuclear thickness. Apparently, the 
A-dependence of the single scattering probability is quite different from A 1 / 3 . In the dipole 
and the BDMPS approach, the accumulated transverse momentum is proportional to the 
mean number of scatterings, i.e. uTa, and hence proportional to A 1 / 3 . Therefore, it was 
concluded in [2] that it is essential to sum all rescatterings in order to get an A 1 /' 3 law. The 
higher twist approach, however, does not only use the double scattering approximation, it 
is also an expansion in gTa- To leading order in this parameter, P\ is identical to the mean 
number of rescatterings. This property of the Poisson distribution is the reason why the two 
expressions for b(j>\) T in Eq. (fTK|) can be so similar. In fact, it has been shown recently ^H] 
that Eq. © remains valid, if the projectile quark exchanges an arbitrary number of gluons 
with the target nucleus. It should be stressed at this point, that 5{p^) n only depends on 
the average color field strength in the target and is not sensitive to details of the color field. 
Regarding details of the pr dependence of nuclear effects, one certainly has to expect very 
different expressions from the dipole approach and the higher twist formalism. 



3 Phenomenological applications 

One can now choose a particular model of the dipole cross section to get an estimate for q-ji 
and Alqs- In this paper, the parametrization of Kopeliovich, Schafer and Tarasov (KST) ^B] 
will be used, because it is motivated from the phenomenology of soft hadronic interactions. 
With the KST-parametrization, C^(0) = C^(0, s) depends on the energy E p of the projectile 
parton, s = 2m^E p , where is the nucleon mass. In all calculation, we also take into 
account Gribov's inelastic corrections (i.e. gluon shadowing), as explained in [21 • At fixed 
target energies, this leads only to a ~ 10% reduction of CV(0, s) for a heavy nucleus, but at 
larger values of y/s, which are relevant for LHC, C^(0,s) is reduced by approximately 1/3. 

Fig. |21 shows the energy dependence of nuclear broadening for quarks and of the three 
parameters CV(0, s), q^ = 9p^Cy(0, s)/2 and 

^'" %,%.(W W ^ 4 (16) 

where (Ta) = J d 2 bT\{b) / A is the nuclear thickness function averaged over impact parameter 
b, and the strong coupling constant is evaluated at a scale S(p'^) :F = (Ta)Cf(0, s). Since this 
scale is in most cases too small for perturbative QCD, we use a running coupling constant 
that freezes at low scales 

Ait 

q1ti / Q 2 +0.54GeV 2 \ 
9 111 0.04GeV 2 ) 

in the spirit of [TJ]. In all three approaches, broadening has a quite significant energy 
dependence, which is due to gluon radiation included in the KST-parametrization. 

We obtain a value of Alq S ~ 0.008 at yjs = 22 GeV (which is the quark energy relevant 
for Fermilab fixed target kinematics) that is very close the one of Ref. jlj (Alq S = 0.01 GeV 2 ), 
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A=200 



dipole/BDMPS 




Figure 2: Upper panel: Broadening 
for a quark propagating through a 
large nucleus as function of y/s = 
y/2m N E p . Because of the hard 
scale in a s , broadening is smaller 
in the higher twist (HT) approach 
(dashed curve) than in the dipole 
and the BDMPS approach (solid 
curve). The dashed curve is cal- 
culated from Eqs. W^TB) with scale 
M = 5 GeV in a s . The other three 
panels show the energy dependence 
of the nonperturbative parameters 
of each approach. 



10 



100 



1000 



V2m N E p (GeV) 



see Fig. El The latter value was determined in Ref. [3] from E772 data on broadening for 
DY. For the gluon transport coefficient, one obtains (^(a/s = 22 GeV) ~ 0.11 GeV 2 / fm, 
more than 2 times as large as the one estimated in Ref. [TJ (q A ~ 0.045 GeV 2 / fm). 

In the dipole approach, broadening only depends on the energy of the parton and not 
on the mass of the dilepton. In the higher twist formalism, however, S{p^) :F depends on the 
dilepton mass through a s . As a consequence, for and Z° production in pA scattering 
with y/S = 8.8 TeV at the LHC {x x « x 2 « 0.01), one has 5(p 2 T )^ ipolc ~ 1.5 GeV 2 for a heavy 
nucleus with A ~ 200 and <5(pt)ht ~ ^-5 GeV 2 . Of course, this estimate assumes that one 
can still apply these formalisms at x 2 = 0.01. As explained in more detail in Ref. jTHj, at 
very low x 2 , the DY cross section is affected by quantum mechanical interferences, and the 
transverse momentum broadening of the produced boson does not reflect the broadening of 
the projectile quark any more. Nuclear broadening in DY at very low x 2 has been calculated 
in Ref. JH] and is expected to be much larger than at medium-low x 2 > 0.01. 

At fixed target energies, however, these interference effects are negligible, and experimen- 
tal data for broadening in DY can be compared to a calculation of broadening for quarks. 
The solid curves in Fig. |3] are obtained from 

HPt)?a ~ S(Pt)?b = ((Ta) - (Tn)) Cfc(0, s), (18) 

where the mean nuclear thickness is calculated with realistic parametrizations of nuclear 
densities from Ref. [25- The dashed curves in Fig. Elare obtained by rescaling o~(Pt)^a by 
the ratio of strong coupling constants, a s (M 2 ) / a s (5(p%) 1 ^ A ). The higher twist formalism is 
strictly speaking not applicable to light nuclei, since all contributions that are not enhanced 
by a power of A 1 / 3 are neglected in this approach. Nevertheless, we believe that a calculation 
with realistic nuclear densities is a reasonable extrapolation to lighter nuclei. 
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Figure 3: Calculations vs. experimental data \2IA \21l HJIij / for broadening with respect to 
(w.r.t.) deuterium (D). Only the N A3 point [231 represents broadening w.r.t. a proton. 
The space in between the curves is an estimate of the theoretical uncertainty. Inner er- 
ror bars show statistical errors, outer errorbars statistical and systematic errors added 
in quadrature. Some points have been slightly displaced for better visibility. The beam 
energies are 800 GeV (Exxx), 140 GeV and 286 GeV (NA10) and 200 GeV (NA3). 



The relevant quark energies for the 800 GeV proton beam at Fermilab are 20 GeV < 
y/s < 25 GeV. The lower solid and dashed curves in Fig. |3] are calculated for y/s = 20 GeV 
and the upper ones for yfs = 25 GeV. For the higher twist calculation, we vary the scale of 
a s in between the J/ip and the T mass. This may serve as an estimate of the theoretical 
uncertainty. As already noted in Ref. [2], the dipole approach overestimates the DY data 
from E772 [20] by several standard deviations. This large discrepancy cannot be explained 
by uncertainties in the parameterization of the dipole cross section [2|. 

However, we point out that the E772 values for broadening [20] were extracted only from 
DY data with transverse momentum p? < 3 GeV [251 > where the pr-differential DY cross 
section is still nuclear enhanced, and may therefore underestimate the true value of S(p'^) :F . 
Moreover, the 0(a s ) parton model does not describe some of the pr-integrated DY cross 
sections measured by E772, either j2Hj- A future analysis j2Z] based on new E866 data |2"H] . 
will include DY data with transverse momentum up to pr < 5 GeV, and may yield values 
of S(p^) :F that are twice as large [2HJ. One can therefore regard the curves in Fig. [21 as 
predictions. 

It is interesting to note that, while E772 only used dileptons with px < 3 GeV, the 
transverse momentum imbalance in photoproduction of dijets was measured by E683 (3U] 
only for jets with pt > 3 GeV. It has been argued in [3T], that the unusually large effect 
observed by E683 is (in part) caused by this restriction on px- In fact, a value of A^qg ~ 
0.1 GeV 2 is needed to accommodate the E683 result 0. The analysis presented in this paper 
clearly favors a much lower value, which is more consistent with the DY data. 

The calculations in the dipole approach for broadening of gluons agree quite well with 
J/i/j and T data, which are underestimated by the higher twist formalism, see Fig. |3] (right). 
Of course, broadening for gluons is equal to broadening in J/ifj and T production, only if final 
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state effects are negligible. This assumption is justified by the observation that broadening 
is very similar (within errorbars) for J/ip and T mesons. 

4 Summary 

In this paper, we quantitatively related the color dipole approach [H |2| to the higher twist 
factorization formalism |3J E] , and studied transverse momentum broadening of fast partons 
propagating through cold nuclear matter. In both approaches, broadening is proportional 
to the average color field strength experienced by the projectile parton We find 

that the corresponding coefficients differ only by the scale of the strong coupling constant. 
While broadening is an entirely soft process in the dipole approach, the extension of the 
QCD factorization theorem to twist-4 is justified by the smallness of a s . In the higher twist 
formalism, a s enters at the typical hard scale of the process that probes the transverse 
momentum of the projectile parton. The equivalence between the dipole and the BDMPS 
approach [Hj was already known before jH]. 

Since the dipole cross section is much better constrained by data than A^qg and q-ji, 
one is now able to obtain new estimates for the latter two quantities. So far, A^qg could 
be determined only from the same data the higher twist approach is supposed to describe 
[I], and estimates for qn were based mostly on physical intuition [XH| - With the KST- 
parameterization of the dipole cross section ^B], which we use, broadening is a function of 
the energy of the projectile parton, as one would expect from a soft process. In the higher 
twist approach, there is an additional scale dependence through a s . To our best knowledge, 
this is the first time that quantitative results for the energy dependence of A^qg and qn are 
presented. It will be necessary to take this energy dependence into account, when applying 
the higher twist formalism and the BDMPS approach at RHIC or even at LHC energies. 

At fixed target energies, numerical calculations in the dipole approach exceed results 
obtained in the higher twist formalism by a factor of ~ 2. Most importantly, the uncer- 
tainty bands of both approaches do not overlap, if one varies the remaining free parameters 
within reasonable limits. Available data, however, do not yet allow to rule out one of the 
theories. Though the dipole approach describes J /if) and T data well, this agreement has 
to be interpreted with great care, since final state effects are not taken into account by the 
theory. We argue, however, that the similarity between broadening for J/ip and T mesons 
indicates that final state effects are rather small. The higher twist approach underestimates 
broadening for J /if) and T mesons. Broadening for DY, on the other hand, is overestimated 
in the dipole approach, while the higher twist formalism reproduces these data well. How- 
ever, the small values of 6(p'^) :F measured by E772 may be the result of a too low p T cut 
imposed on the data. A reevaluation of the E772 data in question, as well as new results 
from E866 measurements, are expected soon |2Zj. This new analysis will probably yield 
significantly larger broadening for DY dileptons [23 I2H]- We stress that no parameter in 
our calculations has been adjusted to fit the data. Thus, the curves presented here can be 
regarded as predictions. 
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